A Universal Magnification Theorem for Higher-Order Caustic Singularities 



A. B. AazamQ 

Department of Mathematics, Duke University, Science Drive, Durham, NC 27708 

A. O. PetterfQ 

Departments of Mathematics and Physics, and Fuqua School of Business, 
Duke University, Science Drive, Durham, NC 27708 

We prove that, independent of the choice of a lens model, the total signed magnification always 
sums to zero for a source anywhere in the four-image region close to swallowtail, elliptic umbilic, and 
hyperbolic umbilic caustics. This is a more global and higher-order analog of the well-known fold 
and cusp magnification relations, in which the total signed magnification in the two-image region 
' of the fold, and the three-image region of the cusp, are both always zero. As an application, we 

, construct a lensing observable for the hyperbolic umbilic magnification relation and compare it with 

■ the corresponding observables for the cusp and fold relations using a singular isothermal ellipsoid 

^SJ ' lens. We demonstrate the greater generality of the hyperbolic umbilic magnification relation by 

showing how it applies to the fold image doublets and cusp image triplets, and extends to image 
configurations that are neither. We show that the results are applicable to the study of substructure 
on galactic scales using observed quadruple images of lensed quasars. The magnification relations are 
also proved for generic 1-parameter families of mappings between planes, extending their potential 
range of applicability beyond lensing. 
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One of the key signatures of gravitational lensing is the occurrence of multiple images of lensed sources. The 
' magnifications of the images in turn are also known to obey certain relations. One of the simplest examples of 
a magnification relation is that due to a single point-mass lens, where the two images of the source have signed 
magnifications that sum to unity: + ^2 = 1 (e.g., Fetters et al. 2001 [1, p. 191]). Witt & Mao 1995 [1] generalized 
this result to a two point-mass lens. They showed that when the source lies inside the caustic curve, a region which 
gives rise to five lensed images, the sum of the signed magnifications of these images is also unity: = 1, where 

' , /ii is the signed magnification of image i. This result holds independently of the lens's configuration (in this case, 
the mass of the point-masses and their positions); it is also true for any source position, so long as the source lies 
00 , inside the caustic (the region that gives rise to the largest number of images). Further examples of magnification 
> relations, involving other families of lens models [N point-masses, elliptical power-law galaxies, etc.), subsequently 
followed in Rhie 1997 [1], Dalai 1998 Q, Witt & Mao 2000 H, Dalai & Rabin 2001 [6^, and Hunter & Evans 2001 0. 
More recently, Werner 2008 Q has shown that the relations for the aforementioned family of lens models are in fact 
topological invariants. 

Although the above relations are "global" in that they involve all the images of a given source, they are not universal 
because the relations depend on the specific class of lens model used. However, it is well-known that for a source 
near a fold or cusp caustic, the resulting images close to the critical curve are close doublets and triplets whose signed 
magnifications always sum to zero (e.g, Blandford & Narayan 1986 [15j], Schneider & Weiss 1992 \\M,, Zakharov 1999 
[3, [E Chap. 9]): 

/^i + = (fold) , 
Ml + A*2 + = (cusp) . 

These magnification relations are "local" and universal. Their locality means that they apply to a subset of the total 
number of images produced, namely, a close doublet for the fold and close triplet for the cusp, which requires the 
source to be near the fold and cusp caustics, respectively. Their universality follows from the fact that the relations 
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hold for a generic family of lens models. The higher-order caustics beyond folds and cusps that we consider are five 
generic caustic surfaces or big caustics occurring in a three-parameter space. Slices of the big caustics give rise to five 
generic caustic metamorphoses (e.g., [ij, Chapters 7 and 9). All five caustic metamorphoses occur in gravitational 
lensing (e.g., Blandford 1990 Fetters 1993 [TT'|, Schneider, Ehlers & Falco 1992 [Tz], and fT]). In addition, the 
magnification relations for folds and cusps have been shown to provide powerful diagnostic tools for detecting dark 
substructure on galactic scales usin g qu adruple lensed images of quasars (e.g., Mao & Schneider 1998 [23|; Keeton, 
Gaudi & Fetters 2003 and 2005 [ll,ll|). 

The aim of this paper is to show that invariants of the following form also hold universally for lensing maps and 
general mappings with higher-order caustic singularities: 

i 

In particular, we show that such invariants occur not only for folds and cusps, but also for lensing maps with elliptic 
umbilic and hyperbolic umbilic caustics, and for general mappings with swallowtail, elliptic umbilic, and hyperbolic 
umbilic caustics. Specifically, we prove that the total signed magnification of a source at any point in the four-image 
region of these higher-order caustic singularities, satisfies: 

Ml + M2 + M3 + fM = ■ 

As an application, we use the hyperbolic umbilic to show how such magnification relations can be used for substructure 
studies of four-image lens galaxies. 

The outline of the paper is as follows. Section HIl reviews the necessary lensing and singular-theoretic terminologies 
and results. Section IIIII states our main theorem, which is for generic lensing maps and general mappings. In 
Section llVl the magnification relations are shown to be relevant to the study of dark substructure in galaxies. We 
also employ a singular isothermal ellipsoid lens to compare the hyperbolic umbilic relations to the fold and cusp ones. 
The proof of the main theorem is quite long and so is placed in Appendices |X] and |B] 

II. BASIC CONCEPTS 
A. Lensing Theory 

We begin by reviewing the necessary lensing and singular-theory terminologies. The spacetime geometry for grav- 
itational lensing is treated as a perturbation of a Friedmann universe by a "weak field" spacetime. To that end, we 
regard a gravitational lens as being localized in a very small portion of the sky. Furthermore, we assume that gravity 
is "weak" , so that near the lens it can be described by a Newtonian potential. We also suppose that the lens is static. 
Respecting these assumptions, the spacetime metric is given by 

aoL = - (l + ^) -'dr^ + (l - ^) + (^^' + Bin^^^ V)) , 

where r is cosmic time, (j) the time-independent Newtonian potential of the perturbation caused by the lens, k is the 
curvature constant, and {R,0,(p) are the coordinates in space. Here terms of order greater than 1/c^ are ignored in 
any calculation involving </). 

The above metric is used to derive the time delay function Ty : L — > R, which for a single lens plane is given by 

Ty(x) = l|x-y|2_^(x) , 

where y — (si, S2) G S is the position of the source on the light source plane S = M^, x = (u,w) G L is the impact 
position of a light ray on the lens plane L C M^, and ip : L — > R is the gravitational lens potential. As its name 
suggests, the time delay function gives the time delay of a lensed light ray emitted from a source in S, relative to the 
arrival time of a light ray emitted from the same source in the absence of lensing. Fermat's principle yields that light 
rays emitted from a source that reach an observer are realized as critical points of the time delay function. In other 
words, a lensed image of a light source at y is a solution x G L of the equation (gradTy)(x) = 0, where the gradient is 
taken with respect to x. When there is no confusion with the mathematical image of a point, we shall follow common 
practice and sometimes call a lensed image simply an image. 

The time delay function also induces a lensing map r] : L — > S, which is defined by 



X I — > vi^) = X — (grad-(/')(x) 
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We call T7(x) = y the lens equation. Note that x € L is a solution of the lens equation if and only if it is a lensed image 
because (gradry)(x) — r]{x) — y. Critical points of the lensing map q are those x S L for which det(Jac?7)(x) = 0. 
Generically, the locus of critical points of the lensing map form curves called critical curves. The value t7(x) of a 
critical point x under t] is called a caustic point. These typically form curves, but could be isolated points. Examples 
of caustics are shown in the third column of Figure [TJ For a generic lensing scenario, the number of lensed images 
of a given source can change (by ±2 for generic crossings) if and only if the source crosses a caustic. The signed 
magnification of a lensed image x G i of a light source at y = ri{x) e is given by 

^^""^ " det(Jacr7)(x) ' 

where we used the fact that dct(Jac?7) = det(HessTy) for single plane lensing. Considering the graph of the time 
delay function, its principal curvatures coincide with the eigenvalues of Hess Ty (x) . In addition, its Gaussian curvature 
at (x, ry(x)) equals det(HessTy)(x). In other words, the magnification of an image x can be expressed as 

^^""^ " Gauss(x,Ty(x)) ' 

where y = t/(x) and Gauss(x, Ty (x)) is the Gaussian curvature of the graph of Ty at the point (x, Ty(x)). Therefore, 
the magnification relations are also geometric invariants involving the Gaussian curvature of the graph of Ty at its 
critical points. Readers are referred to [U, Chap. 6] for a full treatment of these aspects of lensing. 



B. Higher-Order Caustic Singularities 



This section briefly reviews those aspects of the theory of singularities that will be needed for our main theorem. 
The central theorem we shall employ is actually summarized in Table H] below. It is also worth noting that the terms 
"universal" and "generic" will be used often. Formally, a property is called generic or universal if it holds for an open, 
dense subset of mappings in the given space of mappings. Elements of the open, dense subset are then referred to as 
being generic (or universal) . See Chap. 8] for a discussion of genericity. 

We saw in the previous section that the time delay function Ty(x), which can be viewed as a two-parameter family 
of functions with parameter y, gives rise to the lensing map r] : L — > M^. The set of critical points of r] consists of 
all X S L such that det(Jac?7)(x) ~ 0. In this two-dimensional setting, a generic lensing map will have only two types 
of generic critical points: folds and cusps (see [l|. Chap. 8]). The fold critical points map over to caustic arcs that 
abut isolated cusp caustic points; e.g, see the astroid caustic in Figure [TJ 

Now, let Tc_y(x) denote a family of time delay functions parametrized by the source position y and c € K. In 
the context of gravitational lensing, the parameter c may denote external shear, core radius, redshift, or some other 
physical input. The three-parameter family Tc.y(x) gives rise to a one parameter family of lensing maps rj^. Varying 
c causes the caustic curves in the light source plane S to evolve with c. This traces out a caustic surface, called a big 
caustic, in the three-dimensional space M x = {(c, y)}; see Figure [TJ Beyond folds and cusps, these surfaces form 
higher-order caustics that are classified into three universal or generic types for locally stable families f],,, namely, 
swallowtails, elliptic umbilics, and hyperbolic umbilics (e.g., Arnold 1986 Tt'I and fll. Chap. 9]). Generic c-slices of 
these big caustics also yield caustic metamorphoses; see Figure [TJ Note that the point o is a degenerate point of the 
lensing map r]^ on the slice c = 0. 

For the three-parameter family Tc_y(x) of time delay functions, the universal quantitative form of the lensing map 
can be derived locally using ri gid coordinate transformations and Taylor expansions, along appropriate constraint 
equations for the caustics (see [12|, Chap. 6] for details). Table [TJ summarizes the quantitative forms of rj^ for the 
elliptic umbilic and hyperbolic umbilic critical points. The quantitative form for the swallowtail will be dealt with in 
future work. Observe that the elliptic and hyperbolic umbilics for Tc.y (or rj^) do not depend on the lens potential, 
apart perhaps from c in the event that c is a lens parameter. 

One can also consider a general, smooth three-parameter family T'c.s(x) of functions on an open subset of that 
induces a one-parameter family of mappings fc between planes, which are analogs of the lensing map. The universal 
form (also known as the generic or qualitative form) of the one-parameter family fc is obtained basically by using 
differentiable equivalence classes of Tc^s that distinguish c from the coordinates of s, to construct catastrophe manifolds 
that are projected into the space {(c, s)} = MxR^ to obtain local coordinate expressions for fc (e.g., Majthay 1985 [l8[, 
Castrigiano & Hayes 1993 Golubitsky & Guillemin 1973 [13|)- These projections of the catastrophe manifolds are 
called catastrophe maps or Lagrangian maps, and they are differentiably equivalent to fc (see [H, pp. 273-275]). Similar 
to the case for Tc.y and its induced lensing map r]^, a generic family Fc^s and its associated map fc has thre e ty pes 
of caustic singularities beyond folds and cusps: swallowtails, elliptic umbilics, and hyperbolic umbilics (e.g., |l7l |. [l|. 



4 



Type 


Big Caustic 


Caustic Metamorphosis 


Swallowtail 


/l 




X 


Elliptic Umbilic 


C=0 




A 

A ° 


Hyperbolic Umbilic 




+ 


C=0 i( a- 

WW 


1 

> 



FIG. 1: The swallowtail, elliptic umbilic, and hyperbolic umbilic are higher-order caustics shown as surfaces or big caustics in 
the three-parameter space {(c, y)} (middle column). Each c-slice of a big caustic yields caustic curves, which for generic slices 
evolve according to the metamorphoses in the rightmost column. The point o occurs for the slice c — 0. In the case of the 
hyperbolic umbilic. note that the a and b caustic curves are exchanged when c varies through c — 0. This classification is due 
to Arnold 1986 [i3|. 

Chap. 9]). The generic forms of fc about fold, cusp, elliptic umbilic, hyperbolic umbilic, and swallowtail singularities 
are shown in Table HI A detailed treatment of these issues is given in [3, [13, [13 ■ 

In summary, the central result about caustic singularities that we shall use can be stated as follows: 

• A generic, smooth three-parameter family of time delay functions Tc.y can be transformed in a neighborhood 
of a caustic into one of the forms in the second column of Table [J using rigid coordinate transformations that 
distinguish c from the component parameters of y P, . 

• A generic, smooth three-parameter family of general functions i^c,s(x), which need not be a time delay family, 
can be transformed in a neighborhood of a caustic into one of the forms in the third column of Table [T] using 
coordinate transformations distinguishing c from the parameters of s [l|, [l^ . 
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Caustic 


Quantitative Lensing Map 


Generic Map 


Fold (2D) 


Ty{u, u) = - X • y + iamt^ + ^amu^ 

+ 2Q^112W 11 + 2^122™ + ga222l' 

r]{u, v) = (fliiJi + iai22U^ + aii2uv , 

iaii2lt^ + ai22UV + \a222v'^) 


Fb{u, v) — SlU + S2V — iti^ — 

{{u,v) = {u , v^) 


Cusp (2D) 


iy[u,v) — i^y — x-y+2°-ii" + gaiiiit 

+ iaii2lt^U + \ax22Uv'^ + ^02222^^* 
ri{u,v) = (aiilt+ ^0122"^ , 01221^"+ |a2222l'^) 


Fs{u, v) = SlU + S2t; - - fsiu^ - jv* 

f (u, u) = (m , + 11'^) 


Elliptic Umbilic (3D) 


Jc,y(li, j = 2y ^ X ' y + — + 

77^(u, v) — (it^ — , -~2uv + 4cii) 


Fc,b{u, v) — SlU + S2V + c(m^ + V^) 

ic{u,v) = (3«^ — 3ii^ — 2cii , Guv — 2ctj) 


Hyperbolic Umbilic (3D) 


rc,y(u,?;) = iy2 - X ■ y + |(ii^ + i-^) + 2cui; 
J7j.(zi, t;) = (u^ + 2cv , + 2cuj 


Fc,b{u, v) = SlU + 52^ + CUV + U"^ + U'* 

fc{u,v) — (— 3it^ — cv , — 3t;'^ — c?i) 


Swallowtail (3D) 




Fc,s{Uj u) = SlU + S2U — ^S2U*^ — 

-icu^ - 

fc(u, v) = (uu + Cti^ + u'' , u) 



TABLE I: For each type of caustic singularity, the second and third columns show the respective universal local forms of the 
smooth three-parameter family of time delay functions Tc,y and family of general functions Fc,a, along with their one-parameter 
family of lensing maps rj^ and induced general maps fc. For the two-parameter case of the fold and cusp, the constants Oij^ 
denote partial derivatives of rc,y(x) with respect to X = {u,v) = (xi, 0:2), evaluated at the origin: aijk = {d^Tc^y/dxidxjdxk){0). 
These constants do not appear in the quantitative forms of the elliptic and hyperbolic umbilic in the second column. This 
implies that the local behavior of Tc,y and rj^ about an elliptic or hyperbolic umbilic does not depend on the lens potential, 
except possibly through c when c is a lens parameter. We have omitted the quantitative form of the swallowtail for Tc,y and 
rj^ because the proof of its magnification relation will appear in forthcoming work. 



III. MAIN THEOREM 



Consider the universal one-parameter family of lensing maps r]^ in Table [D Let denote a lensed image of a source 
at y, that is, y — ri^{xi), and let /i^ be the magnification of x^, which by H]) is /i^ = 1/ det(Jac ?7^)(xi). For the generic 
mappings fc in Table HI we define the analog of magnification as follows: 



* det(Jacfc)(xi) 

where fc(xi) = s or, equivalently, the point (x^, Fc^si^i)) is a critical point in the graph of Fc.s- 
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Theorem 1. For any of the smooth generic three-parameter family of time delay functions T^^y {or lensing maps 
f]^) and family of general functions F^- s (or general mappings fc) in Table\^ and for any source position y and point 
s in the indicated region, the following results hold: 

1. A2 (Fold) Magnification relations in two-image region: 

Ml + M2 = , Ml + 9Jt2 = . 

2. A'i (Cusp) Magnification relations in three-image region: 

Ml + M2 + Ms - , OTi + 9Jt2 + 9}t3 = . 

3. A4 (Swallowtail) Magnification relation in four-image region: 

OTi + OT2 + OT3 + = . 

4- (Elliptic Umbilic) Magnification relations in four-image region: 

Ml + M2 + M3 + /U4 = , Tli+m2+M3+Tl4 = Q . 

5. (Hyperbolic Umbilic) Magnification relations in four-image region: 

Ml + M2 + M3 + M4 = , SWi + Ma + OT3 + = . 

In the theorem, the ^-magnification (resp., 9Jt-magnification) relations are universal or generic in the sense that 
they hold for an open, dense set of three-parameter families Tc^y (resp., general families Fc.s) in the space of such 
families; see and [J, Chaps. 7,8]. Readers are referred to [J, Chap. 8] for a discussion of universality/genericity. 

The magnification relations in Theorem [1] are also geometric invariants. In fact, we saw in equation Q that each fii 
is a reciprocal of the Gaussian curvature. This is also true of the quantities Tli. To see this, recall that the Gaussian 
curvature at the point (xj, -Fc.s(xi)) in the graph of F^^s is given by 

Gauss(x„f^.,(x,))- dct(HessF.,)(xO 



l + |gradF,,,(x,)|2 

But (xi, Fc.s(xi)) is a critical point of the graph, so gradFc^s(x;i) = 0. A computation also shows that 

det(Jacfc) = det(IIessFc,s) • 

Hence 

an,; = ^ 



Gauss(xj,Fc^s(xj)) 



We use the A, D classification notation of Arnold 1973 [21| in the theorem. This notation highlights a deep 
link between the above singularities and Coxeter-Dynkin diagrams appearing in the theory of simple Lie algebras. 
Theorem [T] is also apparently related to a deep result in singularity theory, namely, the inverse Jacobian Theorem 
and its corollary, the Eulcr-Jacobi formula (see Arnold, Gusein-Zade, & Varchenko 1985 [l^). We are thankful to the 
referee for pointing out this link, which is currently being pursued by the authors. 

As mentioned in the introduction, the fold and cusp magnification relations are known P, [lol [isl . [l^ . but we 
restate them in the theorem for completeness. In addition, note that the magnification relation for the swallowtail is 
established only for the generic form; the quantitative lensing case will be taken up in future work. 

The proof of Theorem [1] is very long. Appendix [Xl gives a detailed proof of the M-magnification relations, while 
Appendix [B] provides a proof of the 971-magnification relations. 
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IV. APPLICATIONS 



Before discussing the applications, we recall that the magnification of a lensed image is the flux Fi of the image 
divided by the flux Fg of the unlensed source (e.g., pp. 82-85]): 

Fi 

where the "+" choice is for even index images (minima and maxima) and the "— " choice is for odd index images 
(saddles). Though Fi is an observable, the source's flux Fs is generally unknown. Consequently, the magnification 
/ii is not directly observable and so magnification sums Mi Eire also not observable. However, we can construct an 
observable by introducing the following quantity: 

where the ± choice is the same as above. This quantity is in terms of the observable image fluxes Fi and image signs, 
which can be determined for real systems [H, [ij . 

Now, aside from their natural theoretical interest, the importance of magnification relations in gravitational lensing 
arises in their applications to detecting dark substructure in galaxies using "anomalous" fiux ratios of multiply imaged 
quasars. The setting consists typically of four images of a quasar lensed by a foreground galaxy. The smooth mass 
density models used for the galaxy lens usually accurately reproduce the number and relative positions of the images, 
but fail to reproduce the image flux ratios. For the case of a cusp, where a close image triplet appears, Mao & 
Schneider 1998 [1^ showed that the cusp ^-magnification relation fails (i.e., deviates from zero) and argued that it 
does so since the smoothness assumption about the galaxy lens breaks down on the scale of the fold image doublet. In 
other words, a violation of the cusp magnification relation in a real lens system implies a violation of smoothness in 
the lens, which in turn invokes the presence of substructure or graininess in the galaxy lens on the scale of the image 
separation. Soon thereafter Metcalf & Madau 2001 [2j] and Chiba 2002 [25] showed that dark matter was a plausible 
candidate for this substructure. 

In 2003 and 2005, Keeton, Gaudi & Fetters [Ulil developed a rigorous theoretical framework showing how the fold 
and cusp /i-magnification relations provide a diagnostic for detecting substructure on galactic scales. Their analysis 
employs the i?-quantity ^ for folds and cusps: 

_ + ^J,2 _ Fi- F2 _ ^J.l + ^J.2 + t^a _ Fi- F2 + F3 

where Fi is the observable flux of image i and image 2 has negative parity. For a source sufficiently close to a fold 
(resp., cusp) caustic, the images will have a close image pair (resp., close image triplet); see the close doublets and 
triplets in Figurema,b,d,e). Theoretically, these images should have vanishing i?foid and i?cusp due to the fold and cusp 
magnification relations and so nontrivial deviations from zero would signal the presence of substructure. In [13, 141], 
it was shown that 5 of the 12 fold-image systems and 3 of the 4 cusp-image ones showed evidence for substructure. 

The study above would look at a multiple-image system and consider subsets of two and three images to analyze 
i?foid and i?cusp, respectively. Such analyses are then "local" when more than three images occur since only two or 
three images are studied at a time. Theorem [T] generalizes the above i?-quantities from folds and cusps to generic 
smooth lens systems that exhibit swallowtail, elliptic umbilic, and hyperbolic umbilic singularities. The i?-quantities 
resulting from these higher-order singularities allow one to consider four images at a time and so are more global than 
the fold and cusp relations in terms of how many images are incorporated. The singularity that is most applicable to 
observed quadruple-images produced by the lensing of quasars is the hyperbolic umbilic (cf. Figure [2]). The associated 
i?-quantity is 

_ Ail + M2 + + fM Fi- F2 + Fs - Fi 



+ lAi2l + Ims] + I/m] Fi + F2 + F3 + F4 ' 

where images 2 and 4 have negative parity. 

We now illustrate the hyperbolic umbilic quantity i?h.u. using a well-known model for a galaxy lens, namely, a 
singular isothermal ellipsoid (SIE) lens. The SIE lens potential and surface mass density are given respectively as 
follows: 

ilj{r, ip) = rF{ip) - -r^cos 2(p, K{r, ip) = 
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where F{(p) and G{(p) satisfy G{(p) = F{lp) + F"{tf) by Poisson's equation, and are given explicitly by 

-^cin 



G(^) = 



^Jl — £C0S2(^ 

i?, 



_i / \/2£cos^ \ ^ -1 / \/2£sin(^ 

cos 09 tan + sm tann 

\ ^/l — e cos / \ ^/l — £ cos 2lp 



where i?oin is the angular Einstein ring radius. The parameter e is related to the axis ratio qhy e ^ (1 — (7^)/(l + g^), 
and should not be confused with the cUipticity e = 1 — q. The cusp at </? = is given by 

Using the Gravlens software by Keeton 2001 [2^, we now solve the SIE lens equation for sources on the positive 
horizontal axis in the four-image region of the light source plane, and compute i?h.u.- Let the SIE have ellipticity 
e = 0.35 and shear 7 = 0.05 oriented along the horizontal axis; both of these values are observationally motivated 
[T3l[T3|. Figure [2Ua,b,c) shows three important image configurations for the SIE: the fold, when the source lies close 
to a fold arc and produces a close pair of images about a critical curve; the cusp, when the source lies close to a cusp 
caustic and produces a close triplet of images about a critical curve; the cross-like configuration of four images, when 
the source sits nearer to the center of the astroid-shaped inner caustic curve. Figure [2l[d,e,f) illustrates how the SIE 
image configurations are similar to those of the hyperbolic umbilic lensing map 77^ given in Table HI See Appendix I A 21 
for more on the hyperbolic umbilic ry^. 

We now look at the behavior of i?foid, -Rcusp, and i?h.u. for an SIE. Table HIl compares -Rfoid and i?h.u. for a source 
approaching a fold arc diagonally from the center of the astroid-shaped inner caustic; see Figure [2l^a) . The fold point 
where the diagonal intersects the fold arc is at 

yfold » (0.14055i?ein, 0.14055i?ein) . 

As the source at y approaches yfoid along the diagonal, the values in Table HIl show that i?foid and i?h.u. each approach 
the ideal value of 0, and that i?h.u. approaches i?foid from above. The reason for this is as follows: From Figure [2l^a) 
we see that there are two pairs of images in a hyperbolic umbilic configuration: the fold image doublet straddling the 
critical curve, and whose two images we denote by di,d2, and the pair consisting of the outer two images, which we 
denote by 01,02. The quantity Rh.u. then becomes 

^ iMdil - \^^d2 \ + ImoiI - IM02I . 

As the source approaches yfoid along the diagonal, Table HIl shows that the quantities \ndi \ — |Md2l and |/XoJ — I/X02I 
stay roughly constant, though the individual magnifications vary. In addition, near the fold, we see that |/UdJ + |Atd2l 
dominates |/XoJ -t- |/io2li causing the denominator of i?h.u. to approach + l/idaL which is the denominator of i?foid- 
This leads to 

_ iMdil - |Md2l ^ IM03I - IM04I > iMdll - |Md2l ^ 

l^dil + l^dsl iMdil + lAid2l ~ l^dil + iMdsl 

The net effect is that i?h.u. approaches i?foid from above (at least for the path along the diagonal). Furthermore, since 
the quantity + 1/^^2 1 diverges, we see that both i?h.u. and i?foid approach the magnification relation value of 0. 

Table IIIII compares i?h.u. with i?cusp for a source approaching a cusp along the horizontal axis from the center of 
the astroid-shaped caustic curve; see Figure [Ub,c). For these values of the ellipticity and shear, we see from Q that 
the two cusps on the horizontal axis are located at 

y±,p«(±0.48i?cin,0) . (5) 

The table shows that as the source approaches y^^sp along the horizontal axis, the quantity i?h.u. approaches -Rfoid 
from below. In other words, Rh.u. is smaller than i?foid- To see why this happens, consider the triplet of sub-images 
in Figure [2Kb), which we denote by ti,t2,t^, and the extra outer image, denote by o. With this notation, 

1^*11 + 1^*2 1 + 1^*3 1 + 1/^0 1 
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SIE for e = 0.35 , 7 = 0.05 Hyperbolic Umbilic rj^ for c = 0.2 



(a) fold (d) fold 




(b) cusp (e) cusp 




FIG. 2: The first column shows fold, cusp, and cross-like configurations due to an SIE with ellipticity e = 0.35 and shear 
7 = 0.05 oriented along the horizontal axis (Panels a,b,c). The second column shows the same configurations due to the 
hyperbolic umbilic lensing map 77^ in Table |T] with parameter value c — 0.2 (Panels d,e,f). In each panel, the sub-figure on 
the left depicts the caustic curves with source position (solid box) in the light source plane, while the sub-figure on the right 
shows the critical curves with image positions (solid boxes) in the lens plane. For the hyperbolic umbilic, image parities have 
been indicated through ± in the given regions. Note that the cross-like configuration shown for the SIE is not a perfect cross, 
which would be the case if the source were centered inside the astroid-shaped inner caustic. Also, for the SIE fold and cusp 
configurations, the source is actually located inside (rather than over) the cusped curve of the astroid. 



As the source approaches y^gp along the horizontal axis, the values in Table UTTl of the cusp relation | — \fJ.t2 1 + 1/^*3 1 
are positive. The inclusion of the outer, negative parity magnification then subtracts from that positive value, 
yielding 

(l/^tl I - lMt2 I + lMt3 I) - IMoI < iMtl I - \^^t2 I + If^ts I , 

which implies that 

-^h.u. — -^cusp • 

Furthermore, Table Hill shows that \fio\ grows fainter faster than the value of the signed magnification of the triplet, 
which yields 

iMti I + lAit2 I + iMta I > iMti I - lMt2 I + lMt3 I > IMoI • 
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Source 




-Rh.u. 


IP"1 1 IP"2 1 


1 A^*-* 1 1 1 1 


IP"1 1 1 IPa2 1 


IPOi 1 1 \t^02 \ 


(^U.lU-rtcin 5 U.lU-ftcin ) 


n 1 /I 

U.i4 






1 01 


o.Ol 


/I Q 
4.O0 


^U.li-ftcin jU.ll-ftcinJ 


U.lo 


n 1 s 


1 99 


1 99 




/I 9S^ 


(0.12i?ein,0.12i?ein) 


0.11 


0.15 


1.22 


1.22 


11.55 


4.21 


(0.13J?cm,0.13J?cin) 


0.08 


0.12 


1.22 


1.22 


15.83 


4.15 


(0.14i?,i,,0.14i?,i,) 


0.02 


0.04 


1.21 


1.23 


65.17 


4.081 


(0.1405i?cin ,0.1405i?cin) 


0.008 


0.015 


1.21 


1.23 


156.80 


4.078 



TABLE II: The quantities Rh.u. and i?foid for an SIE with e — 0.35 and 7 = 0.05 oriented along the horizontal axis. The source 
approaches the fold point yfoid ~ (0.140557?oin, 0.140557?oin) diagonally from the center of the astroid-shaped inner caustic. 
The quantity | — |Atd2 I is the difference in the magnifications of the images in the close doublet, while \^oi \ ~ iMog I is the 
difference for the remaining two outer images; cf. Figure [2ja). 



Source 


^cusp 


Rh.u. 


l^ill + lMt2l + 1^*3 1 


l^til - 1^*2! + l^tsl 


IMoiI 


(0,0) (center) 


0.52 


0.23 


8.49 


4.46 


2.02 


(0.10i?oin,0) 


0.41 


0.22 


9.58 


3.94 


1.49 


(0.15i?oin,0) 


0.36 


0.21 


10.57 


3.76 


1.29 


(0.20i?oin , 0) 


0.30 


0.19 


12.02 


3.61 


1.12 


(0.25i?ein,0) 


0.25 


0.17 


14.20 


3.48 


0.98 


(0.30i?oin , 0) 


0.19 


0.14 


17.71 


3.38 


0.85 


(0.35i?oin , 0) 


0.14 


0.10 


24.10 


3.30 


0.74 


(0.40i?ein , 0) 


0.08 


0.07 


39.02 


3.23 


0.64 


(0.45i?oin , 0) 


0.03 


0.02 


111.5 


3.18 


0.55 



TABLE III: The quantities Rh.u. and JJcusp for an SIE with e = 0.35 and 7 = 0.05 oriented along the horizontal axis. The 
source approaches the cusp point y^sp « (0.487?ein, 0) along the horizontal axis from the center of the astroid-shaped inner 
caustic. The quantity |/itj | — \^t2 I + lA'ta I is the signed magnification sum of the cusp triplet, while |/io| is the magnification of 
the outer image; see Figure [Jl^b). 



In other words, as the source approaches y^gp along the horizontal axis, the contribution of the outer image |/io| 
to the numerator and denominator of Rh.u. becomes negligible. The net effect, at least for the given horizontal axis 
approach, is that i?h.u. and i?cusp converge, with Rh.u. approaching i?cusp from below as they both approach the 
magnification relation value of 0. 

Finally, though Rh.u. can approximate i?foid and i?cusp for fold image doublets and cusp image triplets, resp., the 
hyperbolic umbilic magnification relation has a more global reach in terms of the number of images included. This is 
because Rh.u. also applies directly to image configurations that are neither close doublets nor triplets; e.g., to cross-like 
configurations as in Figure [2l^c). For instance, it was determined in [l^ that to satisfy the relation |i?cusp| < 0.1 at 

99% confidence, the opening angle must be 6* ^ 30°. By opening angle we mean the angle of the polygon spanned by 
the three images in the cusp triplet, measured from the position of the lens galaxy, which in our case, is centered at 
the origin in the lens plane. For the SIE cross-like configuration shown in FigurejU^c), the opening angle is 6 ^ 140°; a 
perfect cross, which would be the case if the source were centered inside the astroid-shaped inner caustic, has 6 = 180°. 
In other words, to satisfy the cusp relation reasonably well, the cusp triplet must be quite tight as, for example, in the 
SIE cusp triplet shown in Figure [H^b) . By contrast, the quantity Rh.u. applies even for values 6 ^ 30°. (In Table Hill 
note how Rh.u. is smaller than i?cusp for source positions closer to the center (0, 0), which yield more cross- like image 
configurations . ) 

A more detailed study of t he prop erties of Rh.u. would require a separate paper and involve a Monte Carlo analysis 
similar to that employed in [ij, |Tj| to study i?cusp and -Rfoid- Such an analysis of Rh.u. would be applicable to the 
currently known 26 four-image lens systems (courtesy of the CASTLES lens sample [27| ) 

V. CONCLUSION 

We showed that magnification invariants hold universally not only for folds and cusps, but also for swallowtails, 
elliptic umbilics, and hyperbolic umbilics. Specifically, for a source anywhere in the four-image region close to each 
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of these caustic singularities, the total signed magnification is identically zero. This result is universal in that it does 
not depend on the class of lens models used, and is thus an extension of the familiar fold and cusp magnification sum 
relations. We proved that these relations hold for generic one-parameter families of lensing maps with the elliptic 
umbilic and hyperbolic umbilic singularities. We also established the relations for generic one-parameter families 
of general mappings, which need not relate to lensing, for the swallowtail, elliptic umbilic, and hyperbolic umbilic 
singularities. We emphasized that these universal sum relations are geometric invariants, because they are sums of 
reciprocals of Gaussian curvatures at critical points. 

The relevance of these higher order magnification invariants to the study of dark substructure in galaxies was shown. 
Using a singular isothermal ellipsoidal model of a galaxy lens, we constructed a lensing observable for the hyperbolic 
umbilic, denoted i?h.u.j and compared it to the well-known fold and cusp analogues, i?toid and i?cusp- These three 
observables approach their magnification relation value of the closer a source gets to a caustic. Significant deviations 
from this value indicate that the lens in question is not smooth, but has some kind of substructure on the scale of the 
image separations. We showed that i?h.u. is a more global quantity than i?toid and i?cusp because i?h.u. considers four 
lensed images simultaneously, while i?foid considers two and i?cusp three. At the same time, we showed that, as the 
source approached a fold arc or cusp point, the quantity i?h.u. approaches i?foid or i?cusp, respectively. More stringent 
conclusions about the properties of i?h.u. await a full Monte Carlo analysis, akin to the one employed recently to 
examine i?foid and i?cusp- 
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APPENDIX A: PROOF OF THE MAIN THEOREM FOR LENSING MAPS 

1. Elliptic Umbilic 




FIG. 3: Caustic curve for an elliptic umbilic lensing map r)^ in equation (|Aip or Table U The number of lensed images for 
sources in their respective regions is indicated. 

The derivation of the quantitative form of the lensing map in the neighborhood of an elliptic umbilic critical point, 
can be found in pL2l . Chap. 6]. The resulting map is 

2 2 

51 = U — V , 

52 = -2uv + Acv . (Al) 

Here y = (si, S2) is the location of a source on the light source plane S, {u, v) the location of a corresponding lensed 
image on the lens plane L, and c is a constant which signifies that the lens mapping under consideration is one in a 
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one-parameter family of lens mappings. The magnification of an image (u, v) is 



(det(Jacs))-i(w,v) = , , • (A2) 



1 



8cu - 4(m2 + tl2) 

A parameter representation of the critical curve is 

u = c(l + cos(f>) , 
w = csm(f> . 

Inserting these into the lens equation (|A1|) gives the caustic curve: 

si ~ 2c^ cos (j) {1 + cos (j)) , 
S2 = 2c^sin(/i(l — COS0) . 

(Note that our notation differs from that of !J3\.) The caustic curve is shown in Figure [H The region inside the closed 
caustic curve constitutes the four-image region. We now show that for all sources inside this region, the total signed 
magnification is identically zero. 

We begin by considering a special case: sources in the four-image region lying on the horizontal axis; that is, with 
S2 = 0. In this case the lens equation (jAip is solvable. The lensed images are (iy^, O) , (2c, zt^/Ac^ ~ ■^i ) ' 'where 
we note that all four of these images are real because < si < 4c^ inside the caustic curve (see Figure [3]). The 
total signed magnification, obtained by inserting each of these four images into (|A2p and summing over, will be zero. 
For the remainder of this section, therefore, we can restrict ourselves to sources (si, S2) inside the caustic curve with 
S2 ^ 0. Note from the second lens equation (jAip that S2 7^ forces the u-coordinate of each lensed image to be 
nonzero: Vi ^ 0. This fact will prove useful below. 

Let (si, S2 7^ 0) denote the position of an arbitrary source lying off the si-axis inside the caustic curve. Let {ui, Vi) 
denote the corresponding lensed images. The total signed magnification /i at (si, S2) is 



1 1 1 1 



Our goal is to show that this sum is in fact identically zero. Let us begin by eliminating u from the lens equation 
()Aip to obtain a (depressed) quartic in v: 

+ {si^'ic^)v^ +2CS2V- ^ = . (A4) 

Knowing that this quartic must factor as 

{v - vi){v - V2){v - V3){v - V4) ^ , (A5) 

we expand (|A5[) and equate its coefficients to those of (|A4p . As a result we obtain four equations involving the w^: 

vi + V2 + V3 + V4 = , (A6) 
V1V2 -I- V1V3 -I- V1V4 -I- V2V3 + V2V4 + V3V4 = si - 4c^ , (A7) 
V1V2V3 -I- V1V2V4 + V1V3V4 + V2V3V4 = -2cs2 , (A8) 

V1V2V3V4 = ■ (A9) 

Next, we use (|A1[) to express each Ui in terms of Wi, bearing in mind that all Vi ^ 0: 

= ^ • (AlO) 

2wj 

Our procedure is to insert (|A10|) into the total magnification (jA3|) , thereby obtaining an expression involving only the 
Vi, n = n{vi), and to then simplify this expression to zero using (|A6p - (jA9p . 

Unfortunately, the equation /i — fJ-{vi), when written over a common denominator, is quite unwieldy. To simplify 
proceedings, we factor the numerator in terms of powers of S2 : 



s 



{vl + + vl + vj) , (All) 
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— 4cs2 + wii^a + I'lWa + V2V3 + viv^ + V2V^ + v^v^ + V2V4 + v^v^ + viv^ + V2V^ + vav^j , (A12) 

A i f / 4 2 I 24, 42, 42, 24, 24, 42, 42, 42, 24, 24, 2 4\ 
[[V1V2 + Vj^V2 + -WiWa + V2V^ + UjUg + V2V^ + Vj^V^ + V2V^ + V^V^ + V^V^ + V2V^ + ) 

'2 / '2 2 22 22 2 

+ 4c [V1V2V3 + V1V2V3 + V1V2V3 + ViV2V4 + V1V2V4 + V1V3V4 + V2V3V4, 

+ ViV^Vi + V2vlv4 + ViV2vl + ViVsvj + V2V3VI)) , (A13) 



ic3/'/'42 , 24 ,4 2, 42,2 4, 24,42 ,24 ,42 ,42 

— IDCS2 ((Wif2f3 + V1V2V3 + V1V2V3 + V1V2V3 + V1V2V3 + V1V2V3 + V1V2V4 + V1V2V4 + + V2V2V4 

+ K1W3W4 + VIV3V4 + vfv2vl + U1W2W4 + vfvsvl + V2V3VI + + V2V3VI 

+ vjv2vj + viV2vj + vjvsvj + vlvsvj + vivlvj + V2vlvj ) 

W4W2W3W4 + Vivlv3V4 + ViV2vlv4 + ViV2V3v'l)) , (A14) 

ic2/'/'442 , 424, 244, 442, 442, 442, 424, 244, 424, 424 
IDS2 ((l'il'2'f^3 + WiW2t'3 + t'it'21^3 + l^ll^2''^4 + ''^l''^3^4 + ^2^3''^4 + ^l''^2''^4 + ''^1^2^4 + '^l''^3''^4 + ''^2'^3^^ 

+ wji;3ti| + f2f3W4) + 4c^ (wf U2f3t'4 + WiW2f3W4 + f 4W2f3W4 + f iW2f3f4 + W1W2V3W4 

+ Wit'|t'|w4 + t'it'2W3^'4 + ^^lW2^3l^4 + Wll^2^^3W4 + ^'l^'2W3^'4 + ^^l^'2^3l^4 ^ , (Al5) 

/442 , 424 ,244 ,44 2,4 42. 442.42 4 

— d4cS2 (W]^ 112^3 "^4 + W4W2f3f4 + f if2't'3'l'4 + l^ll'2''^3^4 + ^^l''^2W3f4 + Wl'y2f3f4 + W4l'2'y3f4 

+ v'lv2V3v\ + v\v2vlv\ + viV2vlvl + vlv2vlvl + vivlvlvl) , (A16) 



64 {vtv'2vH + z;^2'^'3^4' + ^^2«3^4' + ^Mvtvt) ■ (Al7) 

When written over a common denominator, the numerator oi ^ — fJ,{vi) is therefore (jAlip + (|A12p + (IA13P + (|A14p + 
(jA15|) + (|A16|) + (jA17|) . We now proceed to use (jA6p - (jA9p to simphfy each of these terms, beginning with (jAlip . the 
Sj-term. 



The s|-term. We use (M\i and (|T7)) : 

{Vi + V2 + V3+ 114)^ — 2(wiW2 + V1V3 + V1V4 + V2V3 + V2V4 + V3V4) ~ vl +V2 + vl + vl 

= - 2(si - 4c2) = 8c^ - 2si . (A18) 

The S2-term thus simphfies to 

s^(8c2 - 2si) = 8c^s^ - 2sis^ . (A19) 
Because now there is no fi-dependence, this term has been fuUy simphfied. 



The s|-term. We use 



{Vi + -^2 + W3 + l^i) iviV2 + V1V3 + V1V4 + V2V3 + V2V4 + V3V4) - 3 {V1V2V3 + V1V2V4 + V1V3V4 + V2V3V4) 
= vlv2 + V1V2 + vlv3 + vlv3 + Vivl + V2V'1 + v'lv4 + i;2l'4 + vlv4 + i;ii;4 + 7;2f4 + V3v'l 

= (si - 4c2) - 3(-2cs2) = 6cs2 . (A20) 
The S2"term thus simphfies to 

- Acsl{Qcs2) = -24c^s^ . (A21) 



The Sg'term. We proceed in steps. First, 



{vi + U2 + 1'3 + ■(^4) (Wl^'2t^3 + Wl^'21'4 + Wl1'3t^4 + W2t^3^'4) - 4uit;2W3W4 

= V1V2V3 + V1V2V3 + VlV2vl + V1V2V4 + V1V2V4 + V1V3V4 + vlv3V4 + Vivlv4 
+ V2v1v4 + VlV2vl + VlV3vl + W2W3W4 



S 



(-2cs2)-4( ) (A22) 
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Second, 



Third, 



[VlV2 + V1V3 + VlVi + V2V3 + V2Vi + V3V4) — 2 [V1V2V3 + V1V2V3 + V1V2V3 + ViV2V4: + ViV2Vi 
2 2 2 2 2 2 2\ 

22, 22i 22, 22, 22, 22 
— V1V2 + V1V3 + V2V3 + Vj^V^^ + + W3W4 

^ (s, - 4c^)^ - 2{sl) 6 (-|) = (,si - 4cY - I . (A23) 



{V1V2V3 + V1V2V4 + V1V3V4 + W2W3«4)^ " 2 {viV2 + V1V3 + V1V4 + V2V3 + V2V4 + W3W4) (^^1 W2l^3W4 ) 



222, 222, 222i 222 
= Wll'2^'3 +^'l'"2f'4 +^^1^^3^4 +^2'^3«4 



o2 



= i~2cs2f - 2(,si - 4c2) (^-| j = 2c2.^ + ^ . (A24) 
Fourth, wc combine (jAlSp . (IA23p . and (IA24P as follows: 

r 2 I 2, 2, 2\/'22i 22, 22, 22, 22, 2 2\ o('222i 222, 222, 22 2\ 
(Vi + V2+V3 +U4) [V^V^ +V^V3 +V2V3 +V^V4 + V^V^) - 3 (WiW2«3 +'f^l^^2^'4 + ViV^V4^ +V2V3V4) 

42, 24i 42, 42, 24, 24, 42, 42, 42, 24, 24, 24 
= V1V2 + V1V2 + V1V3 + V2VJ + V1V3 + W2W3 + l'iV4 + ^2^4 + W3W4 + W1W4 + W2''^4 + ^3^4 

= (8c2 - 2.1) (^(.1 - 4cr - y) - 3 (2c'4 + ^ 

2 

= 128c^ - 96c^si + 2Ac^sl - 2s? - lOc^sj - ^ • (A25) 
Finally, using (jA25| and (jA22p . the sf-term simplifies to 



4s^ \^128c^ - 96c^si + 2\c^s{ - 2s\ - lOc^s^ - ^ j + ^c^i^s^) 

= 512c®s^ - 384c^sis^ + 96c^s^s| - 8sfs| - 24c^s| - 2sis| . (A26) 



The S2-term. Once again we proceed in steps. First, we note that 

4c^ {V1V2V3V4 + Vivlv3V4 + ViV2vlv4 + V1V2V3VI) = 4c^liiZ)2W3l'4(wi + ^2 + W3 + ^^4) = 



so that the Sj-term reduces to 



16c S2 (l^l't'2''^3 + t'it'2W3 + l'il'2l'3 + tilt'2W3 + Wit'2W3 + 'f^lW2^3 + '^l''^2''^4 + '^l''^2''^4 + ^1^3''^4 + '^2^3''^4 
+ l'iW3W4 + ^'2^3^4 + vfv2vl + Viv'^v'l + v\v3v'l + V2V3VI + W1W3W4 + V2V3VI 
+ vlv2vj + ViV2vj + vfv3vj + V2V3vj + Vivlvj + V2vlvj) . 



Second, we multiply our equation by 1 = ijit^at'sfl ' ^^^'^^ ^ valid operation since each Vi ^ 0, and group together 
terms with a common denominator to obtain 

5 / VS^l + V2V3 + V2vl + V2V4 + V^V4 + V2vl vfv3 + VlV^ + W4 + V^V4 + V1V4 + V3V4 



4CS2 

V ^'1 V2 

_^ V1V4 + V2V4 + vfv2 + VlV^ + vfv4 + V2V4 ^ vfv2 + V1V2 + V4V3 + V^V3 + + V2V3 \ (A27) 

V3 V4 J 

Third, we use (|ll8)) . (|T7| . and (|J22|) to obtain 

[vi + 1^2 + ^3 + ^4) {viV2 + V1V3 + V1V4 + t;2t'3 + y2V4 + V3V4) - {vjv2V3 + t;i'y2'"3 + ViV2vl 
+ V4V2V4 + V1V2V4 + vlv3V4 + V2V3V4 + Vivlv4 + V2vlv4 + W1W2W4 + V1V3VI + V2V3VI) 
= vlv2 + V1V2 + vfv3 + V2V3 + Vivl + V2vl + vfv4 + 172 1^4 + ^31^4 + ViV^ + V2vl + V3VI 

= (8c2 - 2si)(-4c2 + Si) -sl = -2(4c2 - si)^ - sj . (A28) 
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For each of the four terms in (IA27p . we use ljA28p to simplify it. For example, the first term in (|A27p simplifies as 
follows: 

vsvl + v^V3 + + V2V4 + vlv4 + V2vl _ — 2(4c^ — si)^ — si — vfv2 — V1V2 — vfvz — wivf — vfv4 — vivl 

Vl Vl 

= ^ - vfv2 - ^'2 - vfv^ -v^- vfv4 - vl . 

Vl 

Likewise with the remaining terms in (|A27[) . so that the S2-term reduces to 

Acsl ( (-2(4c2 - sif - 4) (— + — + — + —] - ^1^2 - vjvs - vlv4 - vivl - vlvz 

\ \Vl V2 V-i V4J 

— V2V4 — — — — ViV^ — W2W4 — W3«4 — 3 (Ui + ^2 + ^2 + "2) ) • (A29) 

Fourth, we use (jAlSll . (|A6]l . and (|A20| as follows: 

(yl + vl + vl+ vl) {vi +V2+v^ + V4) 

— {vlv2 + Vlvl + vfvs + w|w3 + U1U3 + U2W3 + vlv4 + vlv4 + U3U4 + + V2vl + V^vf) 

= vl + vl + vl + vl = (Sc^ - 2si) - 6cs2 = -6cs2 . (A30) 
Fifth, we use (|X8|) and ((X9|l to obtain 

1111 U2W3t)4 + 'i;it;3W4 + i;iW2W4 + WlW2'y3 8c 



Vl V2 V3 V4 V1V2V3V4 S2 

Finally, we insert (jA30|) and (IA3ip back into (IA29P to obtain the simplified form of the sf-term: 



(A31) 



4c si ( (-2(4c2 - sif -4) {^] - (6CS2) - 3(-6cs2) ) = -1024c^s^ + 512c^sis^ - 64c^s^s^ + 16c^s^ .(A32) 



The Sj-term: Let us begin with the portion of this term with no c-dependence, namely, 

ie 2/442 I 424, 244, 442, 442, 442, 424, 244 
lbS2[Vj^V2V^ + ViV2V^ + ViV2V^ + ViV2V^ + V^V^V^ + V2V^V^ + ViV2V^ + ViV2V^ 

+ vtvlvt + v^vlvt + vlvtvt + vlvlvt) . (A33) 

Analogous to the S2-term above, we begin by multiplying through twice by 1 = i-^v2V3V4 ' then grouping together 
terms with a common denominator, to obtain 

6 f ^2^3 + ^2^i + "1^4 '"l^l + ^1^4 + "^3^4 ^1^4 + ^2^1 + ^1^2 ^1^2 + ^1^3 + ^2^3 
*2 I 712 ' 7?2 ' 712 ' l2 



V^ Vo VS V 
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Next, we use (|A23P on each of the four terms. For example, the first term simplifies as follows: 

+ v^vl + v^vl — — v^v^ — v\v\ + (— 4c^ + si)^ — 53/2 



^1 

(-4c2 + si)2-s2/2 







1 




1 


-1) 






^ 2 

"5 





Likewise with the remaining terms, so that (|A33[) reduces to 

si (-i [vl -^vl + vl + v\) + ({-4c^ + si) 

Next, we use (|A24I) and ((X9|) to obtain 

1111 _ + + + _ 8(4c2 + si) 

9" ~^ 9" ~^ 9" ~^ 9" 9 9 9 9 9 

t;^ 1;^ Kn^zK 4 
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Using this and (jA18|) . we see that (|A33[) simphfies to 

4 (-3(8c^-2.,) + ((-4c^ + .0'-f) {^^^^)) 

= 512c^S2 - 128c^siS2 ~ 32c2s?4 + 8s?s^ - 400^5^ + 2sxs\. (A34) 
There remains the portion of the Sj'tcrm with a factor of c^, namely 

64c S2 (fiW2f3W4 + V-^V2V-iV/^ + Z;]^Z;2f3W4 + l'lW2f3f4 + V-^V2Vi,V/^ + Z;i'i;2l'3l'4 

+ v\viVzv\ + Vxv\vzv\ + WlW2t'3W4 + ^'lW2^'3^'4 + Vxv\vzv\ + WlW2^'3W4) , (A35) 

which factors as 

64c S2 Wlt'2l'3l^4 \ VxV2 + t^lW2 + 173 + ^2^3 + ^X^j, + 1'2W3 + W1W4 + t^2''^4 + 1^3 1^4 + + U2t'4 + V^V^) 

— —16c S2 (WiW2 + V\V2 + W1U3 + ^2^3 + "11^3 + 'V2Vz + V-^V/^ + W2W4 + V-^V^ + Wit;4 + W2W4 + V-iV^) . 

To further simplify this expression, we use (jA18|) . (|A7[) . and (|A22[) : 

+ w| + W3 + W4) (^1^2 + V1V3 + V1V4 + V2V3 + V2V4 + V3V4) — {vlv2V3 + V1V2V3 + ViV2vl + vlv2V4 
+ V1V2V4 + vlv3V4 + V2V3V4 + Vivlv4 + W2W3W4 + ViV2vl + ViV^vl + V2V3VI) 
= vlv2 + V1V2 + vlvs + V2V3 + Vivl + V2vl + vlv4 + V2V4 + w|w4 + Vivl + W2W4 + V^vl 

= (Sc^ - 2si)(si - 4c2) - sl = -32c'* + 160^51 - 2sl - sj . 
Thus (|M5)) simplifies to 

- 16c24(-32c* + 16c2si - 2sl - sj) = 512c^4 ^ 256c^sis^ + 32c^SiS2 + 16c^s^ . (A36) 
Finally, we add (jA34p and (|A36p to obtain the simplified form of the Sj-term: 

1024c^s^ - 384c'*sis^ + 8s^s^ - 24c^s| + 2sis| . (A37) 



The S2-term. First, we factor it as 

V1V2V3 + ViV2V^ + ViV2V^ + V1V2V4 + V1V3V4 + V2V3V4 + V1V2V4 + V1V2V4 
+vlv3vl + V^Vsvl + Vivlvl + V2vlvl) 
= 16c sl {vlv2V3 + vlv2vl + V1V2VI + vlv2V4 + vlvlv4 + V2vlv4 + vlv2vl + WlUj^l 
+vlv3vl + V^Vsvl + Vivlvl + V2vlvl) . 

Second, we multiply our equation by 1 = i-^v2v^v4, ' S^o^p terms together with the same denominator, and use (jA23|) 
to obtain 

4cs2 ( W1W2 + V1V2 + W1W3 + V2V3 + viVj + U2W3 + W1W4 + V2V4 + v^V4 + + i;2W4 + 1^31^4 

9 .9 Sl\ f 1 1 1 1 

(-4c2 + sif + + + 



^Vl V2 V3 V4 

Finally, with the aid of (|A20p and (IA3ip . the S2-term simplifies to 

4c ^6cs2 - f(-4c2 + sif - ("")) = -512c^s^ + 256c*sis| - 32c2s^s| + 400^8^ . (A38) 



The term with no S2-dependence. We factor this term as 



c,( 2222/222, 222, 222, 22 2\ 
04WiW2"3"4 (,«1«2'"3 + «l"2"4 + ^^1«3«4 + '^2'"3'"4j 
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and then use (K9\i and (IA24P to obtain 




) 



(A39) 



We are do ne: for the simpUfied forms of (|XTT|) - (|XT7l) . namely, equations ^MB, dMH), fSMl), (|l32|, (M7\ . (|X38|) . 
and (|A39p . sum to zero: 

(IM) + (1MB + (IM) + (IM) + (1X371) + (IM) + (IM^ 

= (Sc^s^ _ 2sis^) + (-24c2s6) + (512c^s^ - 3840^^314 + OGc^s^s^ - 8s^^4 " '24:c'^4 - 2sisf) 
+ (-10240*^4 + 512c*sis^ - GAc^sjs^ + IGc^s^) + (l024c^4 - 3840^^314 + S44 - "24:0^4 + 2sisf) 
+ (-5120^4 + 2560^^515^ - 32c2s?4 + 400^5^) + (Sc^s^ + 2sis^) = . 
This completes the proof for the quantitative form of the lensing map in the neighborhood of an elliptic umbilic. □ 



The derivation of the quantitative form of the lensing map in the neighborhood of a hyperbolic umbilic critical 
point, can be found in Chap. 6]. The resulting map is 



FIG. 4: Caustic curve for hyperbolic umbilic lensing map r]^ in equation (|A40[l or Table HI The number of lensed images for 
sources in their respective regions is indicated. 

and the corresponding magnification of an image (m, v) is 



2. Hyperbolic Umbilic 



si = + 2cv , 
S2 — f ^ + 2cu , 



(A40) 




(det(Jac s)) ^ {u, v) 



1 



4:{uv — c 



The critical curves are hyperbolas given by 




(A41) 



and the corresponding caustic curve is 




u 



S2 — 2cu H — ^ • 



(A42) 
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The caustic curve is shown in Figured The region "inside the beak" constitutes the four- image region. We now show 
that for aU sources inside this region, the total signed magnification is identicaUy zero. 

First, note from the second lens equation (|A40[) that since Ui{vi) — -^^T^' '^hich is to say, since Vi does not appear 
in the denominator, we do not need to restrict our analysis to the case where Vi ^ 0, as we did with the elliptic umbilic. 
We will therefore let (si, S2) denote an arbitrary source in the four-image region. The total signed magnification fi at 
(si,S2) is 

Kui, V,) = — + — + — ^- + — ^- — ■ (A43) 

A[UiVi — &) \\U2V2 — C) 4(^3113 — c^j 4(1*4114 — c^j 

We begin by eliminating u from the lens equation (jA40p to obtain a (depressed) quartic in v: 

v'^ - 2s2V^ + 8c^v + {si - Ac^si) = . 

Analogous to (|A6[) - (|A9[) and (jB6[) ~ (jB9[) . we obtain four equations involving only the u,;: 

vi + V2 + V3 + V4 = , (A44) 

V1V2 + V1V3 + V1V4 + V2V3 + V2V4 + V3V4 = -2s2 , (A45) 

V1V2V3 + V1V2V4 + V1V3V4 + V2V3V4 = — 8c^ , (A46) 

V1V2V3V4 — -4c^si + S2 . (A47) 

We then insert Ui{vi) into the total magnification (jA43p and factor the numerator of the resulting expression in powers 
of S2: 

c si {viV2V'i + V1V2VA + f 1W3W4 + ^'2^31^4) , (A48) 

— Cs\ (4c^ [viV2 + ViVz + V2V3 + V1V4+ V2V4 + V3V4) + {vfv2V4 + V1V2V4 + W3W4 + V2V3V4 

+ Vivlv4 + V2V^V4 + ViV2vl + ViVsvl + V2V3VI + vfv2V3 + V1V2V3 + ViV2vl)) , (A49) 

CS2 (I2c^ {vi + V2 +V3 + V4) + 4c^ {vlv2 + vivl + vfvs + V2V3 + vivl + V2vl + vlv4 + v^V4 + vlv4 

+ Vivl + V2vl + V3vl) + {vfv2V3 + vlv2vl + ViV^vl + vlv2V4 + W1W3W4 

+ vlvlv4 + vfv2vl + viv^vl + vfv3vl + vlv3vl + vivlvl + V2vlvl)) , (A50) 

- 32c^° - I2c^ {vl + vl +vl+ vf) - 4c'* {vlvl + + vlvl + vlvl + vlvl + v^vl) 

- c {vlvlvl + vlvlvl + vlvlvl + vlvlvl) . (A51) 

The S2-term. Using (jA46|l . this term simplifies to 



cs 



l{~2s2) = -Sc^'sf . (A52) 



The Sj-term. Analogous to (|A18p . we have v\ + V2+ v1+ vl = 4s2, which, when combined with (jA44p . (|A46|) . and 
fCTT]) . yields 



{yl + ^2 + I3 + vl) {V1V2V3 + V1V2V4 + V1V3V4 + V2V3V4) — (^'l^'2^'3^'4)(^'l + V2 + V3 + V4) 

— vlv2V4 + V1V2V4 + vlv3V4 + V2V3V4 + Vivlv4 + V2vlv4 + ViV2vl + ViVsvl + V2V3VI 
+ vlv2V3 + Vivlv3 + ViV2vl = (4S2)(-8C^) - {~4c^Si) = -320^52 . 



The S2"term thus simplifies to 

-csl (4c^ (-2s2) - 32c^S2) = 40cSl . (A53) 
The S2-term. We proceed in steps. First, analogous to (|A22[) . we have 



V1V2V3 + WiW2^3 + ^^11^2^3 + ^'l^'2^'4 + W1W2W4 + WiW3t'4 + ^2^31^4 + t'lW3l'4 
+ V2vlv4 + ViV2vl + V1V3VI + V2V3VI 

= -4(-4c2si -I- sj) = Wc^si - Asl , 
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which, as with (|A28p above, is used to obtain 

(wj + V2+vl+ vj) {viV2 + V1V3 + ViVi + V2V3 + V2V4 + - {vlv2V3 + t'lvf^s + ViV2vl + vIv2Va 

+ Vivjvi + vjvsVi + V^VsVi + Vivjvi + V2vlv4 + ViV2vl + ViVsvj + V2V3vf) 
— vlv2 + V1V2 + vlv3 + V2V3 + Vivl + V2vl + vlv4 + V2V4 + w|w4 + Vivl + V2vl + V^vf 

= (4s2)(-2s2) - (-4s2 + lec^si) = -Asj - IGc^si . (A54) 
Second, analogous to (|A20p and (|A23[) . we have 

V1V2 + V1V2 + + V2V3 + viv^ + V2V3 + v^Vi + V2V4, + v^Vi + viv^^ + V2V^ + V3V^ — -3(— 8c ) = 24c (A55) 

and 

vfv2 + vfvl + v^vj + v\v\ + v\v\ + vlv\ = Asl + 2(-4c^si + si) = &sl - Sc^si . 
Third, we use (|X47l) and ((X55|) to obtain 

viV2VzVi {y'lv2 + V1V2 + vlvz + V2VZ + vivl + V2v1 + + vivi + w|w4 + vivl + V2v'l + v^vl) 

= V'lvlvzV^^ + v\v2V3Vi + 1)^1)2^3^4 + Wiw|w|v4 + v1v2vlv4^ + ViV^V^V/^ 
+ vlv2Vzv'l + Vivlvj,v'l + ViV2vlv1 + Wii;2'y3?^4 + Vivlv'ivl + t)lW2W3W4 

= (S2 ~4c2si)24c^ . 
Fourth, we combine these three equations to obtain 

{vivl + vivl + V2v1 + vivl + 1)2 W4 + vivl) (D1W2W3 + Vll'2W4 + ViV^V4, + V2V3Vi) — {vlvlv3V4 + vfy^VsVi + vlv2vlvi 
+ ViV2vlv4 + vjv2vlv4 + ViV^vlvi + vfv2V3vl + V1V2V3VI + ViV2vlvl + vfv2V3vl + Wid|w3d| + ViV2vlvl) 
= vlv2V3 + vlv2vl + V1V2VI + vlv2V4 + W?w|w4 + wiw|w4 + W?W2w| + + Z)?W3w| + V^Vsvl + + t)2f|w| 

= {6sl - 8c2si)(-8c3) - {sj - 4c2si)24c3 = IGOc^si - 72c^s^ . (A56) 
FinaUy, we use (|A47p . (|A54|) . and (|A56p to simphfy the S2-term: 

CS2 (I2c^(0) + 4c3 {-isj - lec^si) + (leOc^si - 72c^s^)) = 96c^SiS2 - 88c^s| . (A57) 



The term with no S2-dependence. Again, we proceed in steps. First, analogous to (jA24p and (IA30p . we have 
vfy^vl + v\v\v\ + v\vlv\ + v\vlvl = (-Sc^)^ - 2{s\ - sy){~2s2) = 64c^ + 4s2(-4c2si + s\) 

and 

v\^v\ + vl + w| = -24c^ . (A58) 

Second, we have 

{v\v2V3 + Vxv\v3 + VxV2v\ + v\v2V4 + V\v\vi^ + v\v3V4 + v\v3V4 + t)iw|D4 + V2v\v4^ + V\V2v\ + Ulfswl 
+ W2l'3l'4) (^'1^'2 + VYV3 + U2W3 + Ull'4 + V2V4 + 1)31)4) - 3 {v\v\v\ + v\v\v\ + Z)iW3W4 + w|l)3l)4) 
— 3wit)2W3W4 {v\ + 1)2 + 1)3 + 1)4) — 3t)ii;2W3t'4(wiW2 + V\V3 + t)2t)3 + V^V^ + W2f4 + t'3l^4) 
= v\v\v3 + D4w|w3 + v\v2v\ + WiDjfs + t)iW2l'3 + V\v\v\ + 1)^1)2^^4 + v\v\vi4 + t)if|w3t)4 

+ wi^w|i)4 + v\v2v\v/^ + V]v\v\v4 + u| 1)31)4 + ujt;|t;4 + u|u|t)4 + v\v2v\ + Uif|w4 

+ v\v3v\ + v\v2V3v\ + DiW2Z)3W4 + W2l'3W4 + ^'lW2^'3^'I + vxv\v\ + 1)21)31)4 

+ v\v2v\ + t)ll)2l)| + t)iW3^'4 + V2V3v\ + DlWgwl + V2W3WI 

= (-4^2 + 16c2si)(-2s2) - 3 (64c^ + \s2{-^c^ai + si)) - 3(s^ - 4c2si)(4s2) - 3(s^ - 4c2si)(-2s2) 

= -192c^ + 40c2siS2- lOs^ , 
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which we use to obtain 



{vfvl + vfvl + V^vl + vfvl + V^vl + v'^vl) {VIV2 + V1V3 + V2V3 + V1V4 + V2V4 + V3V4) 

— (Wiw|w3 + V1V2V3 + vlv2V^ + ViV^V^ + ViV2vl + ViV^vl + vlv2V4 + vlv2V4 + VIV2V3V4 

+ vlvlv4 + v1v2v'^V4 + ViV2v'^V4 + u|w|ti4 + W^u|f4 + u|u|u4 + Uj*U2l'4 + V1V2V4 

+ vfvsvl + vlv2V3vl + Vivlvsvl + V2V3VI + viV2vlvl + ViV^vl + V2V^vl 

+ v1v2vl + ViV^vl + ViV^vl + V2V3VI + ViV^vl + V2vlvl) 

= vlv2 + + V2V3 + vfvl + v^vl + v^vl 

= (4si + 2 {-Ac^si + sfj) (-2S2) - (-192cS + A0c^siS2 ~ lOs^) 

= 192c6 - 24:C^siS2 ~ 2sl . (A59) 

Third, we have 

{V1V2V3 + V1V2V4 + V1V3V4 + W2W3W4)(wiW2 + V1V3 + V1V4 + V2V3 + V2V4 + V3V4) -~ 3(ui +V2+V3+ U4)(wiW2W3W4) 

= vfvlvs + vlv2vl + vivjvl + vlvlv4 + V1U3U4 + vlvlv4 + v1v2vl + Viv^vj + vlv^vj + V2V3VI + vivjvl + W2W3W4 
= (-8c3)(-2s2) = 1QC^S2 , 

which we use to obtain 

{vfy^vl + vfvlvl + vfvlvl + (wiW2«3 + V1V2V4 + V1V3V4 + V2V3V4) - V1V2V3V4 {y\v\v3 + v\v2v\ + v\v\v\ 

,.2„,2„, I „,2„,2„, I „,2„,2„, I „,2„, „,2 , „, „,2„,2 , „,2„, „,2 , „,2„, „,2 , „, „,2„,2 , „, „.2„,2\ 



+ v\v\v4 + v\v'\v4 + v\v\v4 + v\v2v\ + Vxv\v\ + v\v3v\ + v\v3v\ + Vxv\v\ + l'2w|w|) 

= (64c6 + 4s2 (-4c2si + sD) (-8c3) - (-40^51 + sDiX^e' S2) 



= -512c^ + 192c^siS2 - 480^5^ . (A60) 

Finally, we use (jA58[) . (jA59[) . and (jA60[) to simplify the term with no S2-dependence: 

- 32c^° - 12j(-243) - ic^ (I92c^ - 2Ac^siS2 - 2s^) - c (-512c^ + 192c^siS2 - 48c^s^) 

= -96c®SiS2 + 56c^s^ . (A61) 

We can now verify that the simplified forms of (|A48[) - (|A51|) . namely, equations (|A52p . (jA53|) . (jA57p . and (IA6ip . sum 
to zero: 

((J52I1 + ((J53I) + (ETZl) + (IMTjl = {-8c*sl) + (40cS^) + (96c^siS2 - 88cS^) + (-960*^3132 + 56c*s^) = . 
This completes the proof for the quantitative form of the lensing map in the neighborhood of a hyperbolic umbilic. □ 

APPENDIX B: PROOF OF THE MAIN THEOREM FOR GENERIC MAPPINGS 

1. Elliptic Umbilic 

The derivation of the generic form of a one-parameter family of maps between planes in the neighborhood of an 
elliptic umbilic critical point, can be found in [18, Chap. 8]. The resulting map is 

si = 3v^ — 3u^ — 2cu , 

S2 = Guv — 2cv , (Bl) 
and the corresponding magnification of an image (w, v) is 

(det(Jacs))-(u,.) = ,^2_34^^.^) - (^^^ 

(Note that our notation differs from that of [IBl-) The caustic curve is shown in Figure [51 Although these equations 
are noticeably different from their lensing map analogues IjAip and (jA2|) . it is nonetheless true that the total signed 
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FIG. 5: Caustic curve for the generic elliptic umbilic mapping fc in equation HBlf) or Table |T] The number of lensed images for 
sources in their respective regions is indicated. 

magnification for quadruply-imagcd sources is zero in this the generic case as well. The proof, in fact, is virtually 
identical to the lensing map case in Appendix lA 11 as we now show. 

As before, we begin by considering the special case of a source in the four-image region with S2 — 0. In this case the 

generic lens equation (|Bip is solvable. The lensed images are (— c ± Vc" - 3si ) , 0) , , ±y^^^ , all of which 

are real because — c < si < | inside the caustic curve (see Figure [S]). The total signed magnification, obtained by 
inserting each of these four images into (|B2[) and summing over, will be zero. Once again, therefore, we will restrict 
ourselves to sources (si, S2) inside the caustic curve with S2 ^ 0. Note from the second lens equation ()Bip that S2 ^ 
forces the u-coordinate of each lensed image to be nonzero: Vi ^ 0. 

Let (si, S2 7^ 0) denote the position of an arbitrary source lying off the si-axis inside the caustic curve. Let {ui, Vi) 
once again denote the corresponding lensed images. The total signed magnification /x at (si, S2) is 

""'^ " 4c2-36(u2 + i;2) + 4c2-36(ii2 + v|) + 4c2 - 36(ui + w|) ^ 4c2 - 36(^1 + vl) ' ^^^^ 
We begin by eliminating u from the lens equation (jBip to obtain a (depressed) quartic in v: 

v'~\{s,+c:')v^~lcs2V-'l=0. (B4) 
3 9 36 

Knowing that this quartic must factor as 

(v - vi){v - V2){v - vs.){v - Vi) ^ Q , (B5) 
we expand (jB5[) and equate its coefficients to those of (jB4[) . As a result we obtain four equations involving the 



Vl + V2 + V3 + V4 = , (B6) 

V1V2 + V1V3 + V1V4 + V2V3 + V2V4 + V3V4 — --(si+c^), (B7) 

2 

V1V2V3 + V1V2V4 + V1V3V4 + V2V3V4 = gCS2 , (B8) 

V1V2V3V4 = • (B9) 

Next, we use the second equation in (jBip to express each Ui in terms of w^, bearing in mind that all Vi 7^ 0, 

32 + 2cVi 

Ui{vi)^ — • (BIO) 
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Once again, our procedure will be to insert (|B10[) into the total magnification (jB3p . thereby obtaining an expression 
involving only the Vi, ^ = fJ-{vi). When we do so, and factor according to powers of S2, we obtain 

6 (' 2 I 2 I 2 I 2\ 
- 4 [Vi +V2+V3 +W4) , 

22 2 2 22 2 2 2 2 2\ 

— 4CS2 ('yiW2 + V1V2 + V1V3 + V2V3 + V1V3 + V2V3 + + V2V4 + V3V4 + + + V^V^^) , 

A i in / i 2 I 24, 42, 42, 24, 24, 42, 42, 42, 24, 24, 24\ 
- 4 §2 (9 [V1V2 + V^V2 + U1U3 + ^2^3 + ■'^1^3 + ^2^3 + ^1^4 + ''^2 ''^4 + ^3^4 + ^I'^i + "^2^4 + "^3^4) 

+ 4c (W1W2W3 + V1W2W3 + U1W2W3 + W1W2W4 + WlW2'^4 + Wl'y3W4 + 'y2'^3^^4 

+ 'yiw|w4 + W2U3U4 + WlW2w| + W1V3W4 + W2W3UI)) , 

— I6CS2 (9(W]^'(;2^3 + W1W2W3 + WiW2'^3 + ''^1^^2''^3 + ''^l''^2''^3 + ''^l'^2^3 + ''^l''^2"4 + ''^l''^2''^4 + ''^l''^3''^4 + ^2^ 

+ i;Jw3W4 + l'2"3"4 + vtv2v1 + V1V2VI + 'yi'y3W4 + V2V3VI + ViV^vl + i;2W3W4 

24 24 24 24 24 2 4\ 

+ W1W2W4 + WlW2'y4 + V1V3V4 + V2V3V4 + V1V3V4 + V2V3V4) 

- 4c^ (wiW2W3«4 + 4uiW2W3«4 + 4uiW2l'3l'4 + 4uiW2«3W4)) , 

I 424i 244, 442, 442, 442, 424, 244, 424, 424 
+ ^1^^2^3 + ^l'^2'^3 + V1V2V4 + V1V3V4 + V2V3V^ + V1V2V4 + V1V2V4 + W1W3W4 + V2V3V4 

+ vfv3V4 + U2''^3l'4) + 4c^ (^vfv2V3V4 + VIV2V3V4 + V1V2V3V4 + ViV2vlv4 + vfv2V3V4 

+ ViV2vjv4 + vjv2V3vl + V1V2V3VI + V1V2V3VI + v1v2V3vj + Vivlv3vj + ViV2vlvj)) , 

f 1 A /442 , 424 , 244 , 44 2, 4 42, 442, 42 4 

— 5184 c S2 (^U]^U2''^3''^4 + V1V2V3V4 + V1V2V3V4 + ViV2V3V^ + ViV2V3V^ + V1V2V3V4 + V1V2V3V4 

24 4,4 24, 424,2 44, 24 4\ 

+ Vj^V2V3V4 + Vj^V2V3V4 + V1V2V3V4 + V1V2V3V4 + V1V2V3V4) , 

— 46656 {viV2v'^v1 + vfy'^v^v'^ + vfv'^v^v'^ + V1V2V3VI) . 

Perusal of these expressions shows that aside from constant factors, the form of each of the expressions in parentheses 
is identical to its counterpart in the lensing map case in Appendix lA II This means that the procedures we employed 
above, (jA18|l - (jA39p . carry through without modification. Of course, we expect different final answers, since the 
right-hand sides of (|B7|) - (jB9p differ from those of (jA7p - (|A9|) . With that said, we can forgo the labor and only state 
the final results: 

The s^-term: -fc^s^ _ l^^^e ^ 
The s^-term: fc^s^ ^ 

The 4-term: -f cS* - Sc^s^sj - Sc'slsj - fsfs* + McSgS _ 2^^^6 ^ 

The s3-term: fc'^sj + if cSis| + f c^s^^ „ i^c2^6 ^ 

The s^.term: ~8c^sj ~ ^c^s^s\ - \^s\s\ + \s\s\ + f c^sf + \s^s% , 

The S2-term: f c^^ + f cSi4 + ^c^s\s% - , 

The term with no S2-dependence: —^c^s^ + \s\s\ . 

We can now verify that these terms sum to zero. This completes the proof for the generic form of a one-parameter 
family of maps between planes in the neighborhood of an elliptic umbilic. □ 



-144s2 (9 (vfvlvl 
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2. Hyperbolic Umbilic 



The derivation of the generic form of a one-parameter family of maps between planes in the neighborhood of a 
hyperbolic umbilic critical point, can be found in 18, Chap. 8]. The resulting map is 



si = Sw — cv , 

and the corresponding magnification of an image (m, v) is 

(det(Jac s))^"^ (u, v) : 



1 



-c^ + 36uv 



(Bll) 



(B12) 



The caustic curve is shown in Figure [6l We will show that the total signed magnification for quadruply- imaged sources 
is zero in this the generic case as well. What is more, just as the proof for the generic elliptic umbilic was virtually 
identical to its lensing map analogue in Appendix lA 11 so will be the case here. 

First, note from the second lens equation (jBlip that since Ui{vi) = , which is to say, since Vi does not 

appear in the denominator, we do not need to restrict our analysis to the case where Vi ^ 0. Once again, therefore, 
we will let (si, S2) denote an arbitrary source in the four-image region. The total signed magnification /i at (si, S2) is 



n{ui,Vi) 



1 



1 



1 



1 



— + 36uiVi — + 2>Qu2V2 —(? + 36M3W3 
This time the (depressed) quartic in v is 



-I- 36u4ti4 



3s| + c^si 
27 



= 



and the corresponding four equations involving only the are 

Vl + V2 + V3 + V4 = 

V1V2 + V1V3 + V1V4 + V2V3 + V2V4 + V3V4 = 

V1V2V3 + V1V2V4 + V1V3V4 + V2V3V4 = 

V1V2V3V4 = 



27 ' 

C^Si 



3si 



(B13) 
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(B14) 
(B15) 

(B16) 

(B17) 




FIG. 6: Caustic curve for the generic hyperbolic umbilic mapping fc in equation (IBllf) or TableU The number of lensed images 
for sources in their respective regions is indicated. 



We now insert Ui(vi) into the total signed magnification (jB13|) and factor the numerator of the resulting expression 
in powers of s^'- 



46656 c s\ (vxv^v-i -I- V1V2V4 + V1V3V4 + V2V3V4) , 



(B18) 
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— csl (2592 {viV2 + V1V3 + V2V3 + viVi + V2V4 + V3V4) + 139968 (vlv2Vi + + vfvsVi + V2V3V4 

+ ViV^V^ + V2V^Vi + ViV2vl + V1V3VI + V2V3VI + vfv2V3 + V1V2V3 + ViV2vl)) , (B19) 



c S2 (108 {vi + V2 + V3+ Vi) + 7776 (^1^2 + + vfv^ + V2V3 + Viv^ + wawf + vfvi + V2V4 + v^v^ 
+ vivl + V2vl + V3vl) + 419904 {vfvlvs + u?W2w| + + + vfvlvi 

+ w|w|w4 + vlv2vl + viv^vl + vfvsvl + v^vsvl + vivlvl + V2vlvl)) , (B20) 



- 4 c^" - 324 {vl + vl+vl + vl) - 23328 (wjuf + w^wf + vlv^ + vlvl + v^vl + v^vl) 

- 1259710 c {v^v^vl + vlvlvl + v^vlvl + v^vlvl) . (B21) 

Perusal of these expressions shows that aside from constant factors, the form of each expression in parentheses is 
once again identical to its counterpart in the Icnsing map case in Appendix lA 21 This means that the procedures we 
employed above, (|A52[) - (jA61[) . carry through without modification. Of course, as with the elliptic umbilic, we expect 
different final answers for each term, since the right-hand sides of (|B15I) - (|B17I) differ from those of (|A45p - (IA47p . With 
that said, we once again forgo the labor and only state the final results: 

The s^-term: 1728c'*s^ , 

The s^-term: -8640c''sf , 

The S2-term: 1728c6siS2 + 19008c'^sl , 

The term with no S2-dependence: — 1728c^siS2 — 12096c'*S2 ■ 



We can now verify these terms sum to zero. □ 



3. Swallowtail 



The generic form of a one-parameter family of maps between planes in the neighborhood of a swallowtail critical 
point can be found in Golubitsky & Guillemin 1973 0, p. 176]. The resulting map is 

si = uv + cu^ + u"^ , 

S2 = V , (B22) 
and the corresponding magnification of an image (u, v) is 

(det(Jacs))-^Ki;) ^ ^cu + lu^ + v ' ^^^^^ 

The caustic curve is shown in Figure [T] The "tail" constitutes the four-image region. We will show that the total 
signed magnification for all sources inside this region is identically zero. 

Let (si,S2) denote an arbitrary source in the four-image region, and {ui,Vi) the corresponding lensed images. The 
total signed magnification /i at (si,S2) is 

iiluA = 5-^3 \ 5-^5 1 3_ \ J_ (B24) 

2cui + 4u!f + S2 2cu2 + 4u^ + S2 2cu3 -h 4w| + S2 2cu4 + 4u| S2 

where we have used the second lens equation (|B22p to substitute Vi = S2 into the magnification of each lensed image. 
This time the (depressed) quartic is in u, 

+ cv? + S2U — Si = , 

and the corresponding four equations involving only the Ui are 

ui + U2 + U3 + U4 = , (B25) 

U1U2 + U1U3 + U1U4 + U2U3 + U2U4 + U3U4 = c, (B26) 

U1U2U3 + U1U2U4 + U1U3U4 + U2U3U4 = -S2 , (B27) 

U1U2V3V4 = — si . (B28) 
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S2 




FIG. 7: Caustic curve for the generic swallowtail mapping fc in equation (|B22|) or Table H] The number of lensed images for 
sources in their respective regions is indicated. Note that in order to produce the tail, we must have c < for the given form 
of the swallowtail in HB22[I . 

In what is now becoming a familiar story, we write (|B24[) over a common denominator and then factor the resulting 
numerator in powers of S2'- 

-^sl , 

si (6c (mi +U2+UZ + ua) + 12 {ul + ul + ul + ul)) , 

8c {uiU2 + UiU^ + 'U2W3 + U1U4 + + W3U4) + 16c (^1^2 + U1U2 + + ^2^3 + tiiMg + 

+ ului + ului + ulu4 + Ulul + U2ul + U3M4) + 32 {ulul + uful + U2U3 + ulul + ulul + U3M4)) , 

8c (U1U2U3 + U1U2U4 + U1U3U4 + M2U3li4) + 16c (m^U2U3 + U1U2U3 + U1M2U3 + UlU2Ui + U1M2U4 

+ U\UzUi + M2U3W4 + M1U3U4 + M2M3U4 + UiM2u| + U1M3UI + U2W3M4j + 32c {u\U2U^ 

+ M^U2W3 + M1U2W3 + M^U2U4 + M^M3U4 + U2m|u4 + uJm2u| + U1M2W4 + u\u^u\ 

+ MaWS""! + '«lU3'"4 + U2ulu^^ + 64 [u\ulul + MiU2'^4 + u\ulu\ + M2"3"4) • 

Perusal of these expressions shows that every polynomial in parentheses has already been calculated in the case of 
the hyperbolic umbilic. Once again, we expect different final answers for each term, since the right-hand sides of 
(IB26[) - (IB28[) differ from those of l|A45[) - (|A47[) . With that said, we once again forgo the labor and only state the final 
results: 

The s^-term: 4sf , 

The s^-term: -36si , 

The S2-term: 80^52 + 32csiS2 + 96s| , 

The term with no S2-dependence: — 8c^S2 — 32csiS2 — 64s2 . 



We can now verify that these terms sum to zero. □ 
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